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Purpose of Selected Summary 



• Review some recent activities in Ames’ computational 
methods development area 

• Place these in historical and current context 

• Material selected is to be viewed as core capabilities for 
potential extension to earth science simulations 

• This is intended for open discussion on issues and options 



Historical Examples 



• Grid generation Methods and Tools 

Developed elliptic, hyperbolic, unstructured point-insertion and local 
optimization, patched and overset (Chimera) grid procedures, Cartesian 
grid procedure. 

Resulting codes include: OVERGRID, 3D GRAPE, HYPGEN, DCF, SAGE, 
SURGRID, PRISIM, DELAUNAY, CART3D . 

• CFD Algorithm Advances 

Developed following algorithms and disseminated throughout the US: 

MacCormack explicit algorithm 

Beam-Warming approximate factorization algorithm 

AF1/AF2 algorithm for small disturbance and potential flow 

Pulliam-Chaussee diagonalized ADI algorithm 

Steger-Warming flux-splitting algorithm 

Pseudo compressibility algorithm for viscous incompressible flow 

Linear-reconstruction unstructured scheme 
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Historical Examples 



• Flow Solvers ard Codes 


Small-disturbance, potential and BL: Bailey-Ballhaus, LTRAN2, BL3D 


TAIR/TWING, TOPS 

Euler/Navler-Stokes: 

Incompressible: 

Aero-elastic: 

Rotorcraft: 

Turbomachinery: 


ARC2D/3D, TNS, PNS, TIGER, CNS , OVERFLOW 

INS3D family of codes 

ENSAERO 

TURNS, OVERFLOW.Rtrcrft 
ROTOR, STAGE 


• Post Processing and Scientific Visualization 

PLOT3D, FAST and others 


^ Historical Examples 


• Application Examples - Aerospace 

- Coupled aerodynamic-structured prediction of flutter 

- Full fighter aircraft performance predictions: transonic cruise (F-16), high 
alpha (F-18 with Dryden), STOVL (Harrier), F-15, F-16 .... 

- Full rotorcraft performance prediction (e.g. V-22) 

- Full simulation of the Space Shuttle stack performance (with JSC) 

- 'Tip-to-tair hypersonic aircraft performance prediction 

- Single- and multi-stage turbomachinery performance prediction 

- Simulation/ redesign of the Space Shuttle Main Engine Hot Gas Manifold 

- Space observation system, SOFIA 

- Space transportation vehicles: X-33, X-34, X-38, Shuttle crew escape 
module 

• Application Examples - Ncn-Aerospace 

- First simulation of an Artificial Heart (Penn State) 

- Development of NASA/DeBakey Left Ventricular Assist Device 

- Simulation of naval vehicles (submarine and propeller with the navy) 




Current Tasks Extendable to Earth Simulation 
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1 0 Geometry Definition and Gnd Generation 

1 . t CAD work 

1 .2 Overset grid generation 

1 .3 Topology determination 

1.4 Auto-scripting 

1 .5 Internal flow gridding issues 

2 0 Overset Grid Connectivity 

2.1 Moving grid capability 

2.2 Hole cutting problem 

Current tools need to be worked on 
DCF / Pegsus / X-ray DCF 

2.3 Minimize CPU and wall clock time 

2.4 Develop parallel version compatible with flow solvers 


Current Tasks Extendable to Earth Simulation 


3 0 Solver 

3.1 Time-accurate algorithm 
Integration scheme 
High-accuracy compact schemes 
Convergence acceleration 

3.2 Parallel codes 

MPI, OpenMP, MLP etc. 

Scalability, interoperability and other computer science issues 

4 0 Intelligent Data Management 

4.1 Feature extraction from unsteady data 

4.2 Data compression for communication at remote sites and storage 

4.3 Data base generation using experimental and computational data 
(soft computing tools to cover wide range of operations) 

5 0 Systems analysis tool for entire rocket engine subsystems 



A INS3D - Incompressible Navier-Stokes Solver 
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• Description 

INS3D is a high- fidelity incompressible Navier-Stokes flow solver 
both for steady and unsteady flows. The objective is to provide CFD 
analysis capability on arbitrarily geometry involving complex flow 
phenomena. 

• Overset and rMock Grid Ootu ns 

- Easier grid generation : modularity, better local grid quality 

- Complex moving body problems 

• Flow Solver 

- Includes different algorithm options 

- Parallel version is being developed 

• Current work 

- Rapid, high fidelity procedure for unsteady turbopump system is 
being developed 

POC: C. Klria* kJriaOnaa. naaa.gov ; D. Kwak, kwakOmall.afC.nasa.gov 


EXAMPLE GRID FOR INS3D 
Shuttle Upgrade SSME-rigl 


Shroud 
grid \ 




Unshrouded Impeller Grid : 

6 long blades 
6 medium blades 
12 short blades 




60 Zones 

14.1 Million Grid Points 

Smallest zone : 74,025 
Largest zone : 899,248 


grid for tip 
clearance 


Overset connectivity is 
currently underway by using 
PEGSUS. 
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A Cart3D : SHARP-L1 Reentry Vehicle Design 
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• Enables rapid conceptual design 

- 1 80 Cases delivered in 1 8 days 
using a desktop machine 

- Test matrix 

Mach 0.3 - Mach 5 
Alpha -5deg - +15deg 
beta -ideg - +1deg 
3 Asymmetric bodyflap deflections 




Simulations performed by S. Pandya, Ames 
Research Center 


Scalability of Parallel Cart3D 
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• Parallel Version of Car3D 

- Parallelization through domain 
decomposition and explicit message 
passing. 

- Domain decomposition technique 
based upon space-filling curves 
permits domain decomposition to be 
performed in parallel and on -the- fly 
(at runtime) 

- Parallel speed-ups of -53 on 64 
processors. 

- Combined with robust multigrid to offer 
exceptionally fast convergence 
across the range of Mach numbers. 

• Validation using Citation Twin-engine 

Business Jet 

- 1 .42M cells, Mach 0.84 alpha 1 .81 deg 
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^ Overset Technology for Complex Configuration 
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• .hinvr v Approach 

- Easier grid generation : modularity, better local grid quality 

- Complex moving body problems 

• O', > V'v \ \-.\s -t h.-.-r 

- Third-order upwind differencing 

- Spalart-Allmaras turbulence model 

- Multi-grid and low Mach number preconditioning 

- Performance: upto 12.3 G FLOPS on 256 node Q2K 

• Rapid, high fidelity solutions for complex geometries obtained: 

- Complete commercial high-lift airliner: within 50 days 
Predicted lift within 2% on approach for a complete Boeing 777-200 
landing configuration in 48 days 

- Part-span flap Trap wing: 18 days 

- Transport aircraft flap redesign : 4 days 

POC: S. Rogers, 4-4481, rogersOnas.nasa.gov 
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Overset Technology for Complex Configuration 




• Overset (Chimera) Grid Approach 

NASA Ames Developed CFD Tools 

- OVERFLOW Navier-Stokes Row Sotver 

- Chimera Grid Tools: Pre- and Post 
Processing 

* Enabling flow simulation technology for 
complex configuration and unsteady 
flow Involving bodies in relative motion 


OVERFIOW+CGT : 1998 NASA Software of the 

Year Honorable Mention 

PLOT30 Visualization Software 

1993 NASA Space Act Award 

FAST Visualization Software 

1 995 NASA Software of the Year Award 
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Overset Technology for Complex Configuration 




• Wide Range of Applications 

- Spacecraft ascent and descent 

- Propulsion 

- Aircraft 

- Hydrodynamics 


♦ Current Development 

- Working toward fully automated grid 
generation 

• Steady and Unsteady capabilities 

- Bodies in Motion, 6 DOF 

- Non-equilibrium chemistry 



Turbopump 


yq Data Compression Using Multi resoultion Algorithm 

• DESCRIPTION 

Numerical solution is, in general, only an approximation to the physical 
solution therefore for many purposes truncation of numerical data is 
acceptable. An efficient data compression algorithm for large 3-D data sets 
is desired for storage and transmission of data. 

• ACCOMPLISHMENTS 

A prototype 3-D truncation compression code has been developed and 
demonstrated using CFD data. This code is based on Beam-Warming’s 
supercompact multiwavelets extended from Harten’s interpolatory method. 

• SIGNIFICANCE 

Since most numerical flow solutions are smooth almost everywhere, this 
algorithm should provide efficient CFD data compression. 

Truncation even with high accuracy (e.g. 1x1 O’ 6 error) may allow significant 
compression thus reducing storage and transmission requirements 

• WORK IN PROGRESS 

Develop a 2nd generation wavelets. 
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Discussion 
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• ARC Strength 
CFD+IT 

• Challenges 
Resources 
Technical 

• Areas of Potential Collaboration 

• Other Items 



